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In this communication, the nonlinear bending behaviour of laminated cylin- 
drical panels in the presence of thermal load has been studied using the finite 
element-based simulation package ANSYS. The efficacy of present numerical 
results is compared with published results. Thereafter, the effect of radius-to- 
span ratio (R/a) on the nonlinear static response of isotropic and laminated 
panels under thermal load is studied here. Moreover, the nonlinear bending 
deflection verses x-axis of cylindrical panels at different thermal loads is pre- 


sented here that will serve as benchmark for future research. 


1. Introduction 


Curved panels are widely used in defence, automo- 
bile, aerospace, civil and other engineering appli- 
cations. Thus, the strength and flexural stiff- 
ness of these thin-wall structures under the ther- 
mal environment is very important for engineers and 
researchers (Thang, Nguyen, and Lee). The gener- 
alized differential quadrature numerical method was 
employed for the nonlinear analysis with thermo- 
mechanical loads of moderately thick laminated 
cylindrical panels and considering the variation of 
thickness for conical panels (Naidu and Sinha). 
The linear bending behaviour of skew plates has 
been studied by employing a shear deformable tri- 
angular finite element (Yoosefian, Golmakani, and 
Sadeghian). They studied the effect of thick- 
ness ratio and lamination sequence for skew plates. 
The geometrically nonlinear bending behaviour of 
composite plates with analytical approach such as 
Navier-type and Levy-type techniques has been 
reviewed by (Madrigal, Navarro, and Chaves). 
The nonlinear dynamic responses of curved pan- 
els have been investigated by employing finite ele- 
ment method with first order shear deformation the- 
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ory considering Von-Karman’s geometric nonlin- 
earity. (Loja, Barbosa, and Scares). The buck- 
ling behaviour of isotropic plates under compres- 
sion and shear has been studied by Qiao, P., & Huo, 
X. (Qiao and Huo) Analytically investigated the lin- 
ear buckling of rectangular plates under shear and 
compression using Rayleigh-Ritz Method (Aghdam 
and Falahatgar), Kantorovich Method (Yuan and Jin 
Shufrin, Rabinovitch, and Eisenberger) and numeri- 
cally by using Finite Strip Method (Smith and Srid- 
haran de Vargas Lisboa and Marczak). The lin- 
ear bending results of isotropic plates employing 
energy method has been by Liew (Saadatpour and 
Azhari), and investigated the effect of boundary con- 
ditions. The static analysis of plates has been per- 
formed by using Galerkin Method (Karami, Shah- 
pari, and Malekzadeh Yu and Shen Kumari). Static 
and Dynamic behaviour of flat panels has been 
examined by using finite element method with First- 
Order Shear Deformation Theory (Kumari and Sax- 
ena ZHAO, LIU, and LI). The nonlinear bend- 
ing response of three-dimensional braided lami- 
nated cylindrical panels has been investigated using 
higher-order plate theory with Von-Karman’s Geo- 
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metric Nonlinearity (Moosavi et al.). Nonlinear 
bending behaviour of functionally graded sand- 
wich panels has been studied by using finite ele- 
ment method with higher-order plate theory (Kumar, 
Narendar, and Gopalakrishnan Saadatpour, Azhari, 
and Bradford). It is noticed that nonlinear static 
analysis of curved panel scarce in literature. Hence, 
author has been studied the nonlinear bending 
behaviour of cylindrical panels in the presence of 
thermal load using commercial finite element soft- 
ware ANSYS. 


2. Numerical Procedure 


The finite element method-based simulation soft- 
ware ANSYS is used to investigate the nonlinear 
bending response of cylindrical panels under ther- 
mal load. Eight-node shell element 281 is used here 
to discretize the cylindrical panels. 

The governing equation for the nonlinear bending 
analysis may be written as: 


[EL ]E}> [Exes @]12} [Eas OL} =F} OD 


Here, K,, is the linear, Ky 7; and Ky 2 are nonlin- 
ear stiffness matrices; {d} is the nodal displacement 
vector represents six degrees of freedom (three dis- 
placement u,, V,, wz; and three rotations q;., q,, qz); 
{F } is the nodal load vectors due to thermal load. 

The Newton-Raphson iterative method has been 
used to solve equation (1). 

Residual Force: 


{AF} =[K, +k. (0 )+K AT )|(F}-[F"} 2) 
x,‘ fao}= {ar} @) 


The iteration process is continued until the differ- 
ence between and becomes lesser than error toler- 
ance (e = 0.001) [1]. In this study error criteria are 
considered by: 
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3. Results and Discussion 


The Figure 1. presents non-dimensional cen- 
tre deflection of fully simply supported (SSSS1) 
anti-symmetric angle-ply [-45/45/-45/45] laminated 
cylindrical panel subjected to a uniform thermal 
load (DT) and the present numerical results are val- 
idated available published results by Toan Thang, 
P., Nguyen (Thang, Nguyen, and Lee) and Swamy 
Naidu, N. v., & Sinha (Naidu and Sinha). In the 
present study, full cylindrical panel is used with the 
following geometric and material properties: 

a/b = 1, a/h = 200, R/a = 10 

E, = 138 GPa, Ep = E3 = 8.28 GPa, Giz = Gi3 = 
Go3 = 6.9 GPa, nyg = Ng = Ng = 0.33, 

a, =0.18 10-° /C, ag = 2710-8 /C 


—O— Present study 
—O- Maleki and Tahani (2012) 
—A— Niuda and Sinha (2006) 


Centeral deflection (u, //)) 





0 200 400 600 $00 
Thermal load (AT “C) 


FIGURE 1. Non-dimensional centre deflection of 
fully simply supported (SSSS1) anti-symmetric 
angle-ply laminated cylindrical panel subjected 
to a uniform thermal load (DT ). 


3.1. Boundary Conditions are: 
SSSSI: u, =u, =qy =0; at x =0, a and 

U, =uz=qx =0; at y =0,b 

Figure 1. illustrates the non-dimensional centre 
deflection (u3/h) of the panel under various uniform 
thermal loads (AT = 0 to 800C) of fully simply 
supported (SSSS1) anti-symmetric angle-ply centre. 
The present numerical results are also validated with 
published FEM results and it is seen that percent- 
age of error is 0.0004% and 0.0005% from previous 
studies. 

Next, studied the effect of radius-to-span ratio on 
nonlinear bending behaviour of simply supported 
isotropic cylindrical panel (a/b = 1, a/h = 100). The 
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FIGURE 2. Nonlinear Bending Response (a) and 
NonlinearBending deflection verses x-axis of sim- 
ply (b) R/a = 5 (c) R/a= 8 (d) R/a =12 


following material and boundary conditions are con- 
sidered here: 

E = 210 GPa, n= 0.3, r = 7800 kg/m’, 

a=2~x 10°°/C 

SSSS82: u,; =u, =u,=0; atx =0,a 

and u, =u, =u,=0; at y =0, b 

It is seen that by increasing the radius-to- 
span ratio (R/a) from 5 to 12 increases the non- 
dimensional center deflection (u3/h) from 1.0 to 1.7 
at thermal load AT = 600 C. Moreover, the Non- 
linear Bending Response and Nonlinear Bending 
deflection verses x-axis of simply (b) R/a = 5, (c) 
R/b = 8. 1s presented in Figure. 2 (b-d) for different 
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(a) R/a=5, u3max/h = 1.058 





(b) R/a=8, usmax/h = 1.421 





(c) R/a = 12, usmax/h = 1.702 


FIGURE 3. Contour plots  oflsotropic 
Cylindrical Panel (a/b = 1, a/h = 100) under 
Thermal Load (a) R/a = 5, u3max/h = 1.058 
(b) R/a = 8, u3max/h = 1.421 (c) R/a = 12, 
u3max/h = 1.702 


radius-to-span ratios (R/a = 5, 8, 12) at various ther- 
mal load (AT = 100.8 C, 200.4 C, 300.0 C, 400.8C, 
500.4 C and 600.0 C). As expected by increasing the 
thermal load increases the non-dimensional central 
deflection (u3/h) of cylindrical panel. 

The Contour plots of Isotropic Cylindrical Panel 
(a/b = 1, a/b = 100) under thermal load (AT = 600.0 
C) for different radius-to-span ratios is shown in Fig- 
ure. 3. It is observed that snap-through thermal 
buckling occurs in deep cylindrical panels. 

Now, studied the effect of radius-to-span ratio on 
nonlinear bending behaviour of simply supported 
angle-ply [45/-45/45/-45/45] cylindrical panel (a/b 
= 1, a/b = 100) under thermal load (AT = 50.4 C, 
100.8 C, 200.4 C, 300.0 C). The following material 
and boundary conditions are considered here: 
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FIGURE 4. (a) Nonlinear bending response of simply supported laminated cylindrical panel (a/b = 1, 


a/h = 100) (b) R/a = 5 (c) R/a = 8 (d) R/a= 12 


E,, = 130 GPa, Eog2 = Es3 = 10.3 Gpa, Gy2g = Gy3 
= Go3 = 6.0 Gpa, 

Ny = N43 = Nos = 0.32, r = 1600 kg/m? 

a, =0.2 x 10-° /C, a2 = 40 x 10~°/C 

SSSS2: u, =u, =u,=0; atx =0, a and y =O, b 

It 1s observed that nonlinear bending behavior 
of laminated composite cylindrical panel is qualita- 
tively similar to the isotropic cylindrical panels. 


4. Conclusion 


The nonlinear bending behaviour of isotropic and 
laminated composite cylindrical panels under ther- 
mal load has been studied in this communication 
using the finite element method based commercial 


software ANSYS. To model cylindrical panel 8 node 
281 shell element is used with 441 elements. Firstly, 
validation study was performed with available lit- 
erature in the thermal environment. Thereafter, the 
effect of radius-to-span ratios is studied here on the 
nonlinear static response of cylindrical panels. It is 
noticed that by reducing the deepness of cylindri- 
cal panels reduces the load carrying capacity and 
increases the non-dimensional central deflection of 
these panels. 
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